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Abstract. Astrophysical reconnection takes place in a turbulent medium. The tur- 
bulence in most cases is pre-existing, not caused by the reconnection itself. The model 
of magnetic reconnection in Lazarian & Vishniac (1999) predicts that in the presence 
of turbulence the reconnection becomes fast, i.e. it is independent of resistivity, but de- 
pendent on the level of turbulence. Magnetic reconnection injects energy into plasmas 
through a turbulent outflow from the reconnection region and this outflow can enhance 
the level of turbulence creating bursts of reconnection. Magnetic reconnection in the 
presence of turbulence can accelerate energetic particles through the first order Fermi 
mechanism, as was predicted in Gouveia dal Pino & Lazarian (2005). We discuss new 
numerical results on particle acceleration in turbulent reconnection, compare the accel- 
eration arising from turbulent reconnection to the acceleration of energetic particles in 
turbulent medium. 



1. Reconnection and Numerical Studies 

It is generally believed that a magnetic field embedded in a highly conductive fluid pre- 
serves its topology for all time due to the magnetic fields being frozen-in (Q]; |28l). Al- 
though ionized astrophysical objects, like stars and galactic disks, are almost perfectly 
conducting, they show indications of changes in topology, "magnetic reconnection", 
on dynamical time scales (see H3). Reconnection can be observed directly in the so- 
lar corona (e.g. El), but can also be inferred from the existence of large-scale dynamo 
activity inside stellar interiors (see |29|) . Solar flares and y-ray bursts (see HI; l34h are 
usually associated with magnetic reconnection. A lot of previous work has concentrated 
on showing how reconnection can be rapid in plasmas with very small collisional rates 
(0; HI), which substantially constrains astrophysical applications of the corresponding 
reconnection models. 

A theory of magnetic reconnection is necessary to understand whether reconnec- 
tion is represented correctly in numerical simulations. One should keep in mind that re- 
connection is fast in computer simulations due to high numerical diffusivity. Therefore, 
if there are situations where magnetic fields reconnect slowly, numerical simulations 
do not adequately reproduce the astrophysical reality. This means that if collisionless 
reconnection is indeed the only way to make reconnection fast, then the numerical sim- 
ulations of many astrophysical processes, including those in interstellar media, which 
is collisional at the relevant scales, are in error. At the same time, it is not possible to 
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conclude that reconnection must always be fast on the empirical grounds, as solar flares 
require periods of flux accumulation time, which correspond to slow reconnection. 

To understand the difference between reconnection in astrophysical situations and 
in numerical simulations, one should recall that the dimensionless combination that 
controls the resistive reconnection rate is the Lundquist numbeiQ, defined as S = Ly^^fk, 
where L x is the length of the reconnection layer, Va is the Alfven velocity, and X = 
r\c 2 /4n is Ohmic diffusivity. Because of the huge astrophysical length-scales L x in- 
volved, the astrophysical Lundquist numbers are also huge, e.g. for the ISM they are 
about 10 16 , while present-day MHD simulations correspond to S < 10 4 . As the numer- 
ical efforts scale as h\, where L x is the size of the box, it is feasible neither at present 
nor in the foreseeable future to have simulations with realistically Lundquist numbers. 



2. Reconnection and Turbulence 

While astrophysical fluids show a wide variety of properties in terms of their collision- 
ality, degree of ionization, temperature etc., they share a common property, namely, 
most of the fluids are turbulent. The turbulent state of the fluids arises from large 
Reynolds numbers Re = LV /v, where L is the scale of the flow, V is it velocity and v 
is the viscosity, associated with astrophysical media. Note, that the large magnitude of 
Re is mostly the consequence of the large astrophysical scales L involved as well as the 
fact that (the field-perpendicular) viscosity is constrained by the presence of magnetic 
field. 

Observations of the interstellar medium reveal a Kolmogorov spectrum of electron 
density fluctuations (see |2j; |4J) as well as steeper spectral slopes of supersonic velocity 
fluctuations (see|20, for review). Measurement of the solar wind fluctuations also reveal 
turbulence power spectrum (24). Ubiquitous non-thermal broadening of spectral lines 
as well as measures obtained by other techniques (see yfc 1 1 If) confirm that turbulence is 
present everywhere we test for its existence. As turbulence is known to change many 
processes, in particular the process of diffusion, the natural question is how it affects 
magnetic reconnection. 

To deal with strong, dynamically important magnetic fields Lazarian & Vishniac 
d2~il henceforth LV99) proposed a model of fast reconnection in the presence of sub- 
Alfvenic turbulence. It is important to stress that unlike laboratory controlled settings, 
in astrophysical situations turbulence is preexisting, arising usually from the processes 
different from reconnection itself. In fact, any modeling of astrophysical reconnection 
should account for the fact that magnetic reconnection takes place in the turbulent en- 
vironment and in most cases the turbulence does not arise from magnetic reconnection. 
The analogy here can be as follows: turbulence that is experienced by the plane does 
not arise from the plane motion, but preexist in the atmosphere. 

LV99 identified stochastic wandering of the magnetic field-lines as the most crit- 
ical property of MHD turbulence which permits fast reconnection and obtained an- 
alytical relations between the reconnection rate and the turbulence intensity and the 
turbulence injection scale. 



The magnetic Reynolds number, which is the ratio of the magnetic field decay time to the eddy turnover 
time, is defined using the injection velocity v; as a characteristic speed instead of the Alfven speed Va, 
which is taken in the Lundquist number. 
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Sweet-Parker model 





Turbulent model 



Figure 1. Left panel. Sweet-Parker model versus the model in LV99. Turbulence 
makes the outflow region much wider and independent of resistivity. From (23). 
Right panel. Acceleration of energetic particles as the magnetic field shrinks as the 
result of reconnection. From Lazarian 2005. 



It worth noting that the LV99 model (see Fig. 1 , left panel) is radically different 
from its predecessors which also appealed to the effects of turbulence. For instance, 
unlike (32) and (|13l) Jacobson (1984) the model does not appeal to changes of the mi- 
croscopic properties of the plasma. The nearest progenitor to LV99 was the work of 
( 25l : l26h . who studied the problem numerically in 2D MHD and who suggested that 
magnetic reconnection may be fast due to a number of turbulence effects, e.g. multi- 
ple X points and turbulent EMF. However, ((H; l26h did not realize the key role played 
by magnetic field-line wandering, and did not obtain a quantitative prediction for the 
reconnection rate, as did LV99. 

LV99 revealed a very intimate relation between turbulence and magnetic recon- 
nection. First of all, it shows that reconnection is a necessary ingredient of MHD turbu- 
lence, this is the process that makes the currently accepted picture of MHD turbulence 
(fl2h self-consistent. Moreover, further research in Eyink et al. (@, henceforth ELV11) 
revealed that the expressions of reconnection rate in LV99 can be obtained from the 
basic fluid turbulence concept of Richardson diffusion. 



3. Reconnection and Plasma effects 

For years plasma effects have been considered essential for fast magnetic reconnection. 
On the contrary, LV99 makes use of the MHD approximation. This raises the issue to 
what extend the LV99 model is applicable to astrophysical plasmas. 

An MHD description of plasmas has been revisited recently in ELV11. There 
three characteristic length-scales were considered: the ion gyroradius p,, the ion mean- 
free-path length £ m fp,i, and the scale L of large-scale variation of magnetic and velocity 
fields. Astrophysical plasmas are in many cases "strongly collisional" in the sense that 
tmfpj Pi, which is the case for the interiors of stars and accretion disks. In such cases, 
a fluid description of the plasma is valid. In the "weakly collisional" £ m f p j S> p,,. The 
ratio 

tmfpi A V A ... 

x i — T — ; (1) 

p,- In A c 

follows from the standard formula for the Coulomb collision frequency (see (? ), 
eq.(1.25)). Here A = 4%rik? D is the plasma parameter, or number of particles within 
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the Debye screening sphere. When astrophysical plasmas are very weakly coupled (hot 
and rarefied), then A is large, e.g. of the order of 10 9 or more for the warm component 
of the interstellar medium or solar wind (see Table 1 in ELV11). For such ratio the 
expansion over small ion Larmor radius p, provides "kinetic MHD equations" which 
differ from the standard MHD by having anisotropic pressure tensor. 

Plasmas that are not strongly collisional can be divided into two cases: "collision- 
less" plasmas for which £ m f p j S> L, the largest scales of interest, and "weakly colli- 
sional" plasmas for which L S> @mfp,i- m the latter case the"kinetic MHD" description 
can be further reduced in complexity at scales greater than £ m f p j (see ELV11). This 
reproduces a fully hydrodynamic MHD description at those scales. For instance, the 
warm ionized ISM is "weakly collisional", while post-CME current sheets and the solar 
wind impinging on the magnetosphere are "collisionless." 

Plasmas that are not strongly collisional further divide into two cases: "collision- 
less" plasmas for which £ m f p j » L, the largest scales of interest, and "weakly colli- 
sional" plasmas for which L 3> &mfp,i- m the latter case the"kinetic MHD" description 
can be further reduced in complexity at scales greater than £mfp,i 

(see ELV11). This 

reproduces a fully hydrodynamic MHD description at those scales, with anisotropic 
transport behavior associated to the well-magnetized limit. Among our examples in 
Table 1 above, the warm ionized ISM is "weakly collisional", while post-CME current 
sheets and the solar wind impinging on the magnetosphere are close to being "colli- 
sionless." 

Additional important simplifications occur if the following assumptions are satis- 
fied: turbulent fluctuations are small compared to the mean magnetic field, have length- 
scales parallel to the mean field much larger than perpendicular length-scales, and have 
frequencies low compared to the ion cyclotron frequency. These are standard assump- 
tions of the (fl2h theory of MHD turbulence. They are the basis of the "gyrokinetic 
approximation" (H3; l3lh . At length-scales larger than the Larmor radius p,, another 
reduction takes place. The incompressible shear-Alfven wave modes exhibit dynam- 
ics independent of compressive motions and can be described by the "Reduced MHD" 
(RMHD) equations (see [12; ID). This fact is essential for the LV99 justifying the use of 
the treatment based on an incompressible MHD model. 

Within the LV99 model, the reconnection rate is determined by large scale mag- 
netoc wandering (see Figure 1), while small scale plasma effects may change the local 
reconnection which are irrelevant for the global reconnection at least in the fully ion- 
ized plasma (seel23h. This conclusion is supported by simulations in ( 14) where plasma 
effect were simulated by using anomalous resistivity. We should mention that although 
plasma effects do not change the global reconnection rate, can be important for the 
acceleration of electrons. 

While Hall MHD is a default for many researchers, ELV1 1 showed that the effects 
of the Hall term on the field wandering is negligible on scales larger than p, even if the 
Hall term in the generalized Ohm equation is dominant. 



4. Reconnection and Particle Acceleration 

Magnetic reconnection results in shrinking of magnetic loops and the charged particles 
entrained over magnetic loops get accelerated (see Figure 1, right panel). This process 
was proposed in de Gouveia dal Pino & Lazarian (6, henceforth GL05) (see also Lazar- 
ian 2005, 2006) for the LV99 reconnection and then was adopted for the collisionless 
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reconnection in The physics of the acceleration is the same although GL05 ap- 
pealed to the 3D magnetic bundles (see Figure 1), while (8) considered 2D shrinking 
islands. The latter is the consequence of the constrained 2D geometry and they present 
a strongly degenerate case in 3D. The difference in dimensions affects the acceleration 
efficiency according to (fl5l) . 

GL05 claimed that the acceleration is of the first order Fermi type. This was tested 
recently in (fl6l ). Below we descibe the numerical set up and the results of calculations. 

In order to integrate the test particle trajectories we freeze in time a data cube ob- 
tained from the MHD models of reconnection and turbulence performed in (1141) and 
inject test particles in the domain with random initial positions and directions and with 
an initial thermal distribution. For each particle we solve the relativistic motion equa- 
tion 

— (ymu) =q(E + uxB) , (2) 
dt 

where m, q and u are the particle mass, electric charge and velocity, respectively, E and 

B are the electric and magnetic fields, respectively, J= (l — m 2 /c 2 ) 1 is the Lorentz 
factor, and c is the speed of light. The electric field E is taken from the MHD simula- 
tions E = —v xB + r\J, where v is the plasma velocity, / = V x B is the current density, 
and r\ is the Ohmic resistivity coefficient. We neglect the resistive term above since its 
effect on particle acceleration is negligible (fl5|) . 




Figure 2. Kinetic energy evolution of a group of 10 4 protons in 2D models of re- 
connection with a guide field B Z =Q.Q and 0. 1 (top panels, respectively). In the bottom 
panel a fully 3D model with initial B z =0.0 is presented. The colors show how the 
parallel (red) and perpendicular (blue) components of the particle velocities increase 
with time. The contours correspond to values 0.1 and 0.6 of the maximum number 
of particles for the parallel and perpendicular accelerations, respectively. The energy 
is normalized by the rest proton mass energy. The background magnetized flow with 
multiple current sheet layers is at time 4.0 in Alfven time units in all models. From 
<© 

In Figure |2j we present the time evolution of the kinetic energy of the particles 
which have their parallel and perpendicular (red and blue points, respectively) velocity 
components accelerated for three models of reconnection. The upper left panel shows 
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the energy evolution for a 2D model without the guide field (as in the models studied 
in the previous sections). Initially, the particles pre-accelerate by increasing their per- 
pendicular velocity component only. Later we observe an exponential growth of energy 
mostly due to the acceleration of the parallel component which stops after the energy 
reaches values of 10 3 -10 4 m p (where m p is the proton rest mass energy). Further on, 
particles accelerate their perpendicular component only with smaller linear rate in a log- 
log diagram. In 2.5D case, there is also an initial slow acceleration of the perpendicular 
component followed by the exponential acceleration of the parallel velocity component. 
Due to the presence of a weak guide field, the parallel component accelerates further to 
higher energies at a similar rate as the perpendicular one. This implies that the presence 
of a guide field removes the restriction seen in the 2D model without a guide field and 
allows the particles to increase their parallel velocity components as they travel along 
the guide field, in open loops rather than in confined 2D islands. This result is reassured 
by the 3D model in the bottom panel of Figure 12 where no guide field is necessary as 
the MHD domain in fully three-dimensional. In this case, we clearly see a continu- 
ous increase of both components, which suggests that the particle acceleration behavior 
changes significantly when 3D effects are considered, where open loops replace the 
closed 2D reconnecting islands. 

5. Conclusions 

The results of these studies can be very briefly summarized as follows: 

1 . Advances in the understanding of magnetic reconnection in the MHD regime, 
in particular, turbulent magnetic reconnection in (20) model motivates the studies of 
whether the reconnection in this regime can accelerate energetic particles. 

2. Contracting magnetic loops in magnetic reconnection in 2D, in the MHD 
regime, provides the acceleration analogous to that observed in PIC simulations, which 
proves that the acceleration in reconnection regions is a universal process which is not 
determined by the details of plasma physics. 

3. Acceleration of energetic particles in 2D and 3D shows substantial differences, 
which call for focusing on realistic 3D geometries of reconnection. Our study also 
shows that apart from the first order Fermi acceleration, additional acceleration pro- 
cesses may occur within reconnection sites. 
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